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ABSTRACT 

Molecular oxygen (0 2 ) has been the target of ground-based and space-borne 
searches for decades. Of the thousands of lines of sight surveyed, only those to¬ 
ward Rho Ophiuchus and Orion H 2 Peak 1 have yielded detections of any statistical 
significance. The detection of the 0 2 Nj = 33 —12 and 54-34 lines at 487.249 
GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to 
a short-lived peak in the time-dependent, cold-cloud 0 2 abundance, while the de¬ 
tection of the 0 2 Nj = 33 - 1 2 , 54 - 34 lines, plus the 7 6 - 56 line at 1120.715 GHz, 
toward Orion has been ascribed to time-dependent preshock physical and chemi¬ 
cal evolution and low-velocity (12 km s -1 ) non-dissociative C-type shocks, both 
of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock 
region that is exposed to a FUV field. We report a re-interpretation of the Orion 0 2 
detection based on new G-type shock models that fully incorporate the significant 
effects the presence of even a weak FUV field can have on the preshock gas, shock 
structure and postshock chemistry. In particular, we show that a family of solu¬ 
tions exists, depending on the FUV intensity, that reproduces both the observed 0 2 
intensities and 0 2 line ratios. The solution in closest agreement with the shock pa¬ 
rameters inferred for H 2 Peak 1 from other gas tracers assumes a 23 km s -1 shock 
impacting gas with a preshock density of 8 x 10 4 cm -3 and G 0 = 1, substantially 
different from that inferred for the fully-shielded shock case. As pointed out previ¬ 
ously, the similarity between the LSR velocity of all three 0 2 lines (~ 11 km s” 1 ) 
and recently measured H 2 0 5 32 -4 41 maser emission at 620.701 GHz toward H 2 
Peak 1 suggests that the 0 2 emission arises behind the same shocks responsible for 
the maser emission, though the 0 2 emission is almost certainly more extended than 
the localized high-density maser spots. Since maser emission arises along lines of 
sight of low-velocity gradient, indicating shock motion largely perpendicular to our 
line of sight, we note that this geometry can not only explain the narrow (< 3 km s _1 ) 
observed 0 2 line widths despite their excitation behind a shock, but also why such 
0 2 detections are rare. 

Subject headings: astrochemistry - ISM: abundances - ISM: individual objects 
(Orion) - ISM: molecules - submillimeter: ISM 
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1. INTRODUCTION 


Accounting for the low abundance of molecular oxygen (O 2 ) in dense (n(H 2 ) > 10 3 cm -3 ) 
molecular clouds has long posed a challenge to both observers and theorists. For almost 30 
years, the search for 0 2 was motivated by predictions of gas-phase chemical models that 0 2 
was a major reservoir of elemental oxygen within dense clouds (e.g. Goldsmith & Langer 1978 


Neufeld, Lepp & Melnic k[1995 ). However, the inability of gas-phase models to account for the 
observed abundance of an increasing number of species highlighted the shortcomings of such 
models. This became evident again following the launch of the Submillimeter Wave Astronomy 
Satellite (SWAS) ( |Melnick etal.||2000| ) and Odin ( |Nordh et~aI7||2003[ ) when their largely unsuc¬ 
cessful searches for 0 2 established upper limits to the 0 2 abundance more than 100 times below 
the predictions of these models (e.g., Goldsmith etal. 2000). Modified chemical models, which 
include the effects of dust grains as sites for molecule formation as well as freeze out, were then 
invoked to explain the low 0 2 abundance and weak emission (cf. Bergin et al.|2000; Hollenbach 
etah| |2009| ). With the launch of the Herschel Space Observatory (Pilbratt et al. |2010[ ), atten¬ 
tion turned to using its greater sensitivity to target regions predicted by these updated chemical 
models to possess high columns of warm 0 2 . 

The Open Time Key Program, Herschel Oxygen Project (HOP), was proposed and selected to 
carry out a survey of warm molecular clouds in the following rotational transitions of 0 2 : (N, J) = 


(3,3) -> (1,2) at 487 GHz, (5,4) -> (3,4) at 11A GHz, and (7,6) -> (5,6) at 1121 GHz (see Gold : 
smith etal. 2011| ). Among the sources of particular interest was the Orion Molecular Cloud 
because of the large column densities of warm molecular gas known to exist toward several 
prominent components, including the Hot Core, the Orion Bar, and H 2 Peaks 1 and 2. While no 
convincing 0 2 emission was detected toward the Hot Core (Goldsmith et al. |2011j ) or Orion Bar 
( jMelnick et al.||2012| ), 0 2 emission has been detected toward H 2 Peak 1 ( [Goldsmith et aH20TT 


Chen et al. 2014). This detection raises two questions: (1) What plausible conditions exist at H 2 
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Peak 1 that could produce the observed 0 2 line intensities and line ratios?; and, (2) Are these 
conditions nonetheless sufficiently rare to explain the absence of detectable 0 2 emission toward 
most other sources searched? 


Intense H 2 emission toward Peak 1 is observed in transitions with excitation energies, E u /k, 
ranging from 1015 to 43,000 K that has long been attributed to collisional excitations due to 
shock waves presumed to result from the high-velocity outflow originating ~ 30" to the SE in 
the BN/KL region ( |Gautier etal. 1976[ Beckwith etal. 1978 Rosenthal, Bertoldi & Drapatz 


2000). The Orion Molecular Cloud is also subject to strong far-ultraviolet (FUV) radiation from 


sources such as 9 1 C Ori (cf. Kristensen etal. 2003). In this paper, we examine the role FUV 
radiation plays in altering the chemical and physical conditions in gas capable of reproducing 


the 0 2 observations reported in Chen etaL[ (20l4). In particular, we focus on FUV-illuminated 
gas subject to the passage of non-dissociative C-type shocks as regions especially conducive to 
the production of elevated 0 2 abundances. 

Fully shielded non-dissociative shocks are generally inefficient at producing 0 2 in the post- 
shocked gas. Even when shock velocities are sufficient to sputter material from grain sur¬ 
faces, all but ~ 1% of the gas-phase oxygen not locked in CO is rapidly processed into H 2 0 
within the warm (T > 400 K) postshock gas via a set of neutral-neutral chemical reactions 
(0+H 2 —* OH+H; OH+H 2 —> H 2 0+H). In the absence of a FUV field, the H 2 0 abundance is 
predicted to remain high throughout the postshock region and the 0 2 abundance never exceeds 
a few percent of that of H 2 0 (cf. praine etal.|1983| Kaufman & Neufeld|| 1996b ). However, the 
presence of a FUV field can significantly increase the postshock 0 2 abundance by: (1) increas¬ 
ing the preshock atomic oxygen abundance through both the photodissociation of gas-phase 
O-bearing species and the photodesorption (from grains) and subsequent photodissociation of 
desorbed O-bearing species, leading to a higher postshock H 2 0 abundance; and, (2) increas¬ 
ing the postshock abundance of OH and O through the photodissociation of postshock H 2 0, 
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enabling increased production of 0 2 via the reaction O+OH —> 0 2 +H. 

In Chen et al. (2014), the three components necessary to model the observed 0 2 emission - i.e., 
the preshock gas, the shock and its postshock thermal profile, and the postshock chemistry - 
were addressed separately. In particular, the preshock gas and shock were assumed to be fully 
shielded, i.e., no FUV radiation reaches these regions, while the postshock chemistry was mod¬ 
eled assuming various levels of FUV illumination. It is unlikely that the preshock, shock, and 
postshock gas, which are in close physical proximity, would experience significantly different 
FUV fluxes; however, the absence of a shock code that incorporates the full effects of a FUV 
field dictated the approach described in Chen et al. 


In this paper, we re-examine the results presented in Chen etal., with several inportant differ¬ 
ences. The main contribution of the present work is the self-consistent treatment of the preshock, 
shock, and postshock regions under the influence of FUV fields common to all components. 
FUV radiation can not only increase the preshock atomic oxygen abundance, as described above, 
but the effects of FUV radiation on the structure of C'-typc shocks can be significant. These ef¬ 
fects include increased postshock temperatures at a given shock velocity, reduced velocities at 
which these shocks break down (and H 2 is dissociated), as well as the altered abundance of key 
postshock molecular species. We will show that: (1) the thermal and chemical changes induced 
by the FUV radiation on the preshock, shock, and postshock gas, and including the sputtering 
of H 2 0 from grain mantles at shock velocities greater than ~ 20 km s -1 , yield multiple combi¬ 
nations of shock velocity and FUV field intensity capable of producing detectable 0 2 emission; 
and, (2) the 0 2 line intensity ratios can be an important discriminator between shock models. 
Use of the observed 0 2 line ratios was not used as a test of the shock models presented in Chen 
etal. 


In addition, the inclusion of a FUV field in all gas components mitigates a problem identified 
in Chen etal., i.e., the need in their model to assume > 10% sputtering efficiency of H 2 0 from 









Page 7 


grain mantles behind a 12 km s -1 C -type shock. As will be discussed, a sputtering efficiency of 
10% or greater behind a 12 km s -1 shock is higher than theoretical predictions as well as recent 
observations. 

Finally, the association between co-located H 2 0 maser emission and the 0 2 emission noted 
in Chen etal. provides important clues regarding the orientation of the H 2 0- and likely the 
0 2 -producing shocks. We point out that this orientation not only justifies the assumption of 
relatively large column densities of 0 2 along the line of sight, as discussed in Chen etal., but 
also provides an alternate explanation to that offerred in Chen et al. for why such 0 2 detections 
are rare. 

In §2, we briefly review the 0 2 observations upon which our analysis is based. In §3, we describe 
the modifications to the 6'-type shock models required when such shocks are exposed to FUV 
radiation. In §4, we describe the assumed shock geometry, the methods used to compute the 0 2 
line integrated intensities, and how well the results fit the observations. In §5, we discuss the 
ways in which FUV-illuminated G'-typc shocks and the shock geometry toward Orion H 2 Peak 1 
combine to provide an explanation for the detected 0 2 emission while also accounting for the 
apparent rarity of such emission. 


2. OBSERVATIONS AND RESULTS 


The Herschel observational results considered in this paper are presented in Chen etal. (2014). 


The interested reader is referred to Chen et al. for details regarding the observations and data 
reduction; only the results we model are summarized here. 


The observations were conducted in 2012 using the Heterodyne Instrument for the Far Infrared 


(HIFI; de Graauw etal. 2010) and were centered at J2000 coordinates R.A. = 5 h 35 m 1452 and 
Dec. = —5° 22' 31". HIFI was used in dual beam switch (DBS) mode with the reference po- 
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sitions located 3' on either side of the source. For each transition, eight local oscillator (LO) 
settings were used to allow sideband deconvolution. The integration time for each LO setting 
was 824 seconds for the 487 GHz and 774 GHz spectra, and 3477 seconds for the 1121 GHz 
spectrum. The observational results relevant to our modeling effort are presented in Table 1. 


3. MODIFIED SHOCK MODELS 

3.1. FUV-llluminated C-Type Shocks 

Within interstellar shocks possessing low fractional ionization and a strong magnetic field, the 
flow variables may remain continuous; i.e., the neutral and ionized gases do not experience dis¬ 
continuities and the gas remains relatively cool and molecular. The physics of G'-typc shocks in 
well-shielded gas has been studied extensively (e.g., Draine|1983 ; Draine et~aL|1983| Kaufman 
& Neufeld 1996a Flower & Pineau des Florets|[2010); however, it is increasingly evident that 


FUV radiation from nearby stars plays a role in the physics of molecular outflows (e.g., Kris 


tensen etal. 2012, van Kempen etal. 2010). In order to assess the role of FUV radiation, we 


have made several modifications to the Kaufman & Neufeld ( |1996a| ) model. These modifica¬ 
tions affect the abundances of important species in the preshock gas, the shock length scale, and 
the abundances of oxygen-bearing species in the postshock gas. We briefly describe the effects 
below; a more detailed description will be included in a forthcoming paper. 


3.2. P re shock Conditions 


To determine the conditions in the preshock molecular gas, we use the detailed photodissoci¬ 


ation region (PDR) model of Hollenbach etal. (2009). This model computes the abundances 


of numerous atomic and molecular species, as well as those of charged and neutral dust grains 
and PAHs, as a function of extinction into a cloud, given the input parameters of gas density, 
7i, and FUV radiation field strength, G 0 (scaling factor in multiples of the average Habing local 
























Page 9 


interstellar radiation field; Habing |1968 ). In order to model a shock propagating through such 
gas, the important outputs of the PDR model are: (i) the types and abundances of charge car¬ 
riers, which determine how well the ionized species couple to the magnetic field; and (ii) the 
abundances of gas-phase species, especially oxygen- and carbon-bearing species that become 
incorporated into important gas coolants. It is the interplay between momentum transfer (moder¬ 
ated by the coupling length) and the efficiency of gas cooling (which limits the gas temperature) 
that determines the final shock structure. As shown by Hollenbach et al. (2009), significant frac¬ 
tions of oxygen and carbon nuclei can be locked up on the surfaces of dust grains in molecular 
gas exposed to moderate FUV fields at extinctions A v ~ 1. Based on the same PDR model de¬ 
scribed in detail in Hollenbach etal. (2009), Fig. [T| shows the preshock gas-phase abundances of 
H 2 0 and O as functions of density and G 0 . 


To first order, the length scale depends only on the ionization fraction and the strength of the 
magnetic field. In gas fully-shielded from FUV radiation, the fractional ionization is set by the 
cosmic ray ionization rate and is low, resulting in a relatively long ion-neutral coupling length 
scale and (7-shocks that can be supported up to velocities of ~ 40 km s -1 . In unshielded gas, the 
fractional ionization reaches its maximum value of ~ 10 -4 , set by the carbon abundance, and 
results in the shortest possible length scale. In gas at Ay = 1, the length scale is between these 
two extremes and depends sensitively on the gas density and grain properties. An important 
point arising from these calculations is that, because the length scale is shorter than in fully- 
shielded gas, the maximum 6'-shock velocity is generally less than 40 km s _1 . 


3.3. Postshock Chemistry 

Within the shock, the abundances of gas-phase O-bearing species are set by competition between 
neutral-neutral reactions, which tend to drive O nuclei into 0 2 , OH and H 2 0, and photodissoci¬ 
ation, which destroys these molecules. Studies of (7-type shocks in fully-shielded gas show that 
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O-nuclei are efficiently driven into H 2 0 once the gas temperature exceeds ~ 400 K, at shock 
velocities >10-15 km s -1 ; in the absence of photodissociation, the H 2 0 abundance remains 
high in the postshock gas. The presence of FUV radiation has two effects. First, it shortens 
the length scale and thus lowers the shock velocity at which the gas reaches 400 K. Second, it 
destroys molecules in the cooling postshock gas; to account for this, we have added photodis¬ 
sociation reactions to the shock chemical scheme, including dissociation of H 2 0, OH and 0 2 . 
Under these conditions, the gas has higher abundances of O and OH than without photodisso¬ 
ciation. The O and OH are driven into molecular oxygen by the reaction O + OH —> 0 2 + H, 
which has a reaction rate that is within a factor of two of 3 x 10“ 11 cm 3 s” 1 over the range 10 K 
to 500 K (McElroy et al .|2013 ). This is in contrast with the neutral-neutral reactions that drive O 
into OH and H 2 0, each of which has barriers of > 1000 K. We find that the postshock 0 2 abun¬ 
dance is directly correlated with the preshock gas-phase O-nuclei abundance; more O-nuclei in 
the preshock gas lead to more H 2 0 in the warm gas followed by more 0 2 in the cooler gas. 

In addition to increases in the preshock gas-phase O-nuclei abundance caused by FUV radiation, 
H 2 0 can be removed from ice mantles by the shocks themselves. Ice mantles on grains are 
sputtered by ion-neutral collisions in shocked regions where there is substantial relative velocity 


between charged grains and neutral molecules. We follow the treatment of Flower & Pineau 
des Florets|(1994) for the removal of oxygen and water-ice mantles, but with a threshold energy 
for ice removal consistent with that found by jNeufeld et al. (2014), i.e. complete grain mantle 
removal for shocks with v , > 25 km s _1 . Note that these results are similar to those found 


by both praine et al.| ( |1983| ) and |Jimenez-Serra et al.| ( |2Q08[ ), who find ~ 10% of mantles are 
removed in 20 km s” 1 shocks. At velocities less than ~ 20 km s _1 , sputtering of grain mantles is 
not expected to be efficient, so the gas-phase abundances of O and C nuclei should be unchanged 
by the passage of such slower shocks. 


Because the 0 2 itself is subject to photodissociation, the 0 2 abundance, N(0 2 )/N(H 2 ) = y (0 2 ), 
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does not reach the abundance of O-nuclei in the gas. An example may be seen Fig. [2] in which 
models with G 0 = 0.1, 1 and G a = 10 are compared. Note that the peak O 2 abundance is higher 
in the G a = 10 case since the higher G 0 keeps more O in the gas phase. Moreover, the higher 
postshock gas temperatures when G 0 is high allow the neutral-neutral reactions to more than 
compensate for the photodestruction of H 2 0 while T > 300 K, keeping the H 2 0 abundance high; 
however, when T drops below ~ 300 K, the H 2 0 photodestruction rate exceeds its formation 
rate and the H 2 0 abundance quickly drops, as shown in Fig. [2] 

The conditions necessary to provide a significant 0 2 column are: (i) a large fraction of O-nuclei 
in the preshock gas; (ii) sufficient FUV radiation to dissociate the shock-produced H 2 0, but not 
so much that no significant 0 2 survives in the postshock gas; and (iii) a FUV flux that remains 
below the threshold at which the length scale will no longer support C'-shocks. 


4. MODELING APPROACH 


4.1. Shock Geometry 


Chen etal. (2014) have identified Orion H 2 Peak 1 as the most likely source of the reported 
0 2 emission. Peak 1 is the site where the high-velocity gas from the dynamical center of the 
outflow, about 30" to the southeast in the vicinity of the BN/KL region, impacts the surround¬ 


ing cloud, creating strong and extended collisionally-excited near-infrared H 2 emission (Cun¬ 


ningham 2006) as well as 22 GHz H 2 0 maser emission (cf. Genzei etal.] 1981). Using Her- 
5c/re//HIFI, Neufeld et al. ( 2013| ) recently reported the detection of H 2 0 maser emission arising 
from the 5 32 - 4 4 i rotational transition at 620.701 GHz toward H 2 Peak 1 (see Fig. [3]). As noted 
in Chen et al., the similarity between the LSR velocities of the maser emission of ~ 11 km s -1 
and that of the 0 2 suggests a connection. 


Along with 22 GHz H 2 0 maser observations obtained during the same epoch, Neufeld etal. in- 
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terpret the maser emission as being consistent with a pumping model in which the population 
inversions arise from the combined effects of collisional excitation and spontaneous radiative 
decay. Moreover, because collisional excitations that produce a population inversion in the 532 - 


4 41 transition require H 2 densities greater than 10 s cm 3 (cf. Neufeld & Melnick 1991), and 
because the upper 5 32 state lies 698 K above the ground state, gas compressed and heated be¬ 
hind a shock front provides the most plausible explanation for the source of this maser emission. 


Finally, the direction of maximum gain for stimulated emission lies along the path with the 
minimum velocity gradient which, for a 67-type shock, is perpendicular to the direction of prop¬ 
agation. Thus, it’s likely that the shocks giving rise to the 22 GHz and 621 GHz H 2 0 masers 
toward H 2 Peak 1 are propagating perpendicular to our line of sight. While we are not sug¬ 
gesting that the 0 2 emission necessarily arises from the H 2 0 maser spots, the presence of the 
H 2 0 masers points to shocks moving in the plane of the sky, some of which achieve high enough 
postshock densities and temperatures to produce maser action, but most of which may not. Such 
a scenario has significant advantages for explaining the 0 2 emission observed toward H 2 Peak 1: 
(1) shocks viewed in cross section (vs. face-on) can potentially present the observer with a large 
column of gas with a relatively high postshock 0 2 abundance; and, (2) the propagation velocity 
of the shocks in the plane of the sky can assume any value (up to the 67-type shock breakdown 
velocity), while the line-of-sight velocity gradient can be low, and thus consistent with the nar¬ 
row (< 3 km s -1 ) 0 2 line widths observed (see Table 1). This geometry is shown schematically 
in Fig. [4} 


4.2. Shock Models and O 2 Emission 

A series of shock models were generated for preshock densities ranging between 10 3 cm~ 3 and 
10 6 cm~ 3 in steps of 0.1 dex and incident FUV fluxes, 67 0 , of 0.1, 1, and 10. For each value of 
G 0 and preshock density, shock velocities ranging from 5 km s _1 up to the 67-shock breakdown 
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velocity, in steps of 1 km s _1 , were also generated. In total, 839 shock models were computed 
for a G 0 = 0.1, 618 models were computed for a G 0 = 1, and 455 models were computed for a 
G 0 = 10. That fewer models were generated as G c increases reflects the fact that as G 0 increases, 
G-type shocks break down at lower shock velocities. Thus, a smaller range of shock velocities 
- and models - are consistent with G-type shocks as G 0 rises. 

For each shock model, the postshock H 2 density, temperature, and H 2 0, OH, and 0 2 abundance 
are computed as a function of distance behind the shock front. The range of postshock distances 
considered extends beyond where: (1) the H 2 0, OH, and 0 2 abundances peak and subsequently 
drop by at least a factor of 10; and, (2) the postshock temperature drops to 10 K. Meeting these 
criteria ensures that our modeling captures essentially all of the 0 2 emission since, beyond 
the postshock distances we consider, the reduced 0 2 abundances combined with the low gas 
temperatures produce negligible contributions to the total 0 2 integrated intensity. 


Within the postshock region in which the density, temperature and 0 2 abundance favor non- 
negligible 0 2 emission (as defined above), the gas is divided into 1000 equally-spaced zones, 
each characterized by an H 2 density, temperature, and 0 2 abundance (relative to H 2 ) as com¬ 
puted by the G-shock code. To obtain the 0 2 line intensities from each zone, the equilibrium 
level populations have been calculated using an escape probability method. We use the rate 
coefficients for collisions between 0 2 and He computed by Lique ( |2010| ), multiplied by 1.37 to 
account for the different reduced mass when H 2 is the collision partner. Since these rates were 
computed for gas temperatures < 350 K, extrapolation to higher temperatures was necessary for 
application to postshock conditions. To do this, a polynomial was fitted to the 13 computed rates 
between 5 K and 350K for each transition and extended to 1400K. As presented in §4.3 (and fig¬ 
ures therein), 0 2 is formed in the postshock gas downstream of where H 2 0 and OH are formed. 
At the postshock distance where the 0 2 abundance approaches its peak, gas temperatures have 
typically dropped below about 500 K. Thus, rates up to 1400 K suffice. 
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Because 0 2 has no dipole moment and can only emit quadrupole radiation, the line center op¬ 
tical depth is small under almost all circumstances. As a result, the emergent line intensity is 
proportional to the 0 2 column density, N( 0 2 ). The range of H 2 densities and 0 2 abundances 
over which the postshock 0 2 emission peaks is determined from the postshock densities, tem¬ 
perature, and abundances computed for each shock model, combined with the radiative transfer 
calculations; however, N( 0 2 ) also depends on the line-of-sight depth of the shock, which is 


difficult to assess. A study of water maser emission from behind C-type shocks by Kaufman & 


Neufeld (1996b) suggests that the aspect ratio of the zone of H 2 0 emission - i.e., the ratio of 


the line-of-sight shock depth to the cross-sectional width of the high-H 2 0 abundance zone on 
the plane of sky - could exceed a factor of 100. For our purposes, we assume the shock width 
(in the plane of the sky) equals the distance between where the postshock gas temperature first 
rises above 10 K to where it first drops back to 10 K. We then multiply this distance by 100 
to obtain the line-of-sight depth of the 0 2 emitting gas. However, our analysis is not tied to 
the H 2 and 0 2 column densities thus derived since the total integrated intensity measured for 
each 0 2 line depends not only on A r (0 2 ), but also on the number of shocked regions within the 
Herschel !HIFI beams, and hence the total spatial extent of the 0 2 emission within these beams. 


The spatial extent of the 0 2 emission was only weakly constrained by Chen etal. (2014), other 
than inferring that the 0 2 emission most likely did not fill the Herschel/lUFl beams at 487, 774, 
and 1121 GHz. Thus, we treat the spatial extent of the 0 2 emission as a free parameter. 


We seek to match both the absolute 0 2 line fluxes and the observed line ratios. We first derive 
the column densities and emission area needed to reproduce the measured 0 2 487 GHz line flux. 
The ratio of the accompanying line fluxes at 774 GHz and 1121 GHz to that at 487 GHz are then 
compared with those observed. We assume that the H 2 and 0 2 column densities and emission 
area needed to reproduce the observed 487 GHz 0 2 line flux are given by: 
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N(R 2 ) 

= 100 x n(H 2 ) dio7 

(1) 

N( 0 2 ) 

= N (H 2 ) x(0 2 ) , 

(2) 


where 100 is the assumed aspect ratio, n(H 2 ) is the H 2 postshock density in the zones where 
487 GHz emission peaks, d 10 is the postshock distance between where the gas temperature first 
rises above 10 K and where it first returns to 10 K, and y(0 2 ) is the 0 2 abundance relative to 
H 2 . Because the 0 2 column density required to reproduce the observed optically thin line flux 
is inversely proportional to the area of the emitting region, we introduce the scaling factor 7 , 
which is defined as 


7 = 


400 sq. arcsec 

^487 


(3) 


where 400 sq. arcsec is chosen to be roughly consistent with the area of the 0 2 emitting region 


assumed by Chen etal. (2014), and 2)437 is the actual area of the 0 2 487 GHz emission in 
sq. arcsec were it known. Since we compute the 487, 774, and 1121 GHz line emission as a 
function of postshock depth, the area of the 774 GHz and 1121 GHz emitting regions, relative 
to 2 f 48 7 , is computed for each shock model and scale as 


Relative area of 774 GHz emitting region 


Relative area of 1121 GHz emitting region 


A 


487 


*487 


w 774 \ 

(4) 

W 4 87 ) 


^1121 \ 

(5) 

1(7487 ) 



where u’ 487 , u'774, uqm are the widths of the 487 GHz, 774 GHz, and 1121 GHz emitting regions 
behind the shock front, respectively. The width of the emitting region in each line is taken to be 
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the range of postshock depths over which the line flux remains greater than 1 percent of its peak 
value. 

Finally, we define a ‘Goodness of Line Ratio Fit’ parameter to measure how closely each model 
comes to reproducing the observed 774 GHz and 1121 GHz line fluxes once the 487 GHz line 
flux has been matched: 


Goodness of Line Ratio Fit 


-R774m — -^7740 

1 

Rll21m ~ Rll21o 

R774o 

1 

Rll21o 


( 6 ) 


where R~ 74 m is the model-predicted ratio of the 774 GHz integrated intensity to the 487 GHz 
integrated intensity, and R 7740 is the observed ratio of the 774 GHz integrated intensity to the 
487 GHz integrated intensity. Likewise, i?n 2 im is the model-predicted ratio of the 1121 GHz 
integrated intensity to the 487 GHz integrated intensity, and i?u 2 i 0 is the observed ratio of the 


1121 GHz integrated intensity to the 487 GHz integrated intensity. Chen etal. (2014) report an 
observed 487 : 774: 1121 GHz integrated line intensity ratio of 1:1.90: 0.53 toward Orion H 2 
Peak 1. As indicated by Eqn. (6), model-predicted line ratios that exactly match the observed 
ratios would result in a ‘Goodness of Line Ratio Fit’ of zero. 


4.3. Shock Models Results 

The results for shock models of varying preshock density, shock velocity and FUV field, com¬ 
puted for G 0 = 0.1, 1, and 10, are shown in Figs. [5] [6j [7J respectively. Only those models that 
produce a ‘Goodness of Line Ratio Fit’ less than 6 for G 0 = 0.1 and 1, and less than 14 for 
G 0 = 10, are shown for clarity. The best-fit model for each value of G 0 is summarized in Ta¬ 
ble 2; the 0 2 line ratios produced by these best-fit models are shown in Fig. [8j The profiles of 
postshock 0 2 line integrated intensity, abundance, and temperature for the best-fit G' (l = 0.1, 1, 
and 10 shocks models are shown in Figs. [9j [T0| and [IT] respectively. Finally, the line center 
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optical depth, t, in the 487, 774, and 1121 GHz transitions was evaluated for each of the best-fit 
results given in Table 2. In no case did r exceed 0.07 and, in most cases, was significantly lower, 
thus justifying the optically thin assumptions made here. 

5. DISCUSSION 


Following more than six years of observations, and more than 7000 Galactic lines of sight sur¬ 
veyed, the SWAS mission reported only one tentative detection of O 2 , toward the Rho Ophiuchi 
molecular cloud (Goldsmith et al .|2002 ). Likewise, the Odin mission reported only upper limits 
to the 0 2 line strengths ( |Pagani etal.||2003[ |Sandqvist etal.||2008| ), with the exception of one 
possible detection, also toward the Rho Ophiuchi cloud (Larsson et a L||2007j ). The upper limits 
to the 0 2 abundance set by SWAS and Odin are more than 100 times lower than that predicted by 
equilibrium gas-phase chemical models. This discrepancy was later understood to be primarily 
the result of the exclusion of dust grains from these models, which serve as important sites for 
both the freeze out of H 2 0 and the surface formation of H 2 0. This water ice subsequently re¬ 
mains locked on grain surfaces until either photodesorbed by FUV or X-ray photons, sublimated 


at grain temperatures above ~ 100K, or sputtered by shocks with velocities > 25 km s' 1 (Draine 
etal. |1983[ Neufeld etal. 2014). By sequestering large amounts of elemental oxygen in water 
Bergin et al. ( ]2000| ) and Hollenbach et al] (}2009j) showed that the gas-phase production of 


ice, 


0 2 is effectively suppressed. 

The Herschel Oxygen Project, guided by these updated models, used HerscheV s greater sensi¬ 


tivity to continue the search for 0 2 . To date, however, Rho Ophiuchi (Liseau etal. 2012) and 
Orion H 2 Peak 1 remain the only sources with statistically significant O 2 detections. It might be 
expected that processes common to both sources, yet rare overall, could account for the 0 2 emis¬ 
sion. Unfortunately, the conditions inferred by Liseau et al. to explain the O 2 emission toward 
Rho Ophiuchi are insufficient to account for the emission toward Orion H 2 Peak 1. In particular, 
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toward Rho Ophiuchi, the 0 2 emission is attributed to a combination of two emitting regions, 
one with iV(0 2 ) > 6 x 10 15 cm' 3 and T < 30 K, and the other with iV(0 2 ) = 5.5 x 10 15 cm -3 
and T > 50 K. The inferred beam-averaged 0 2 abundance is ~ 5 x 10 8 in the warmer com¬ 
ponent, and somewhat higher in the colder component. The successful detection of 0 2 toward 
this source, among the many sources toward which no 0 2 emission was detected, was attributed 
to time-dependent quiescent cloud chemistry - i.e., Rho Ophiuchi was surmised to have been 
observed during a relatively short period when the evolving 0 2 abundance was near its peak. 


Such a scenario is unlikely to apply to Orion H 2 Peak 1. First, the constraints on the spatial ex¬ 


tent of the 0 2 emission toward Peak 1 provided in Chen et al. (2014) - i.e., an 0 2 emitting region 
less than 25" in diameter - require that the 0 2 column density be between about 3 x 10 17 cm' 2 
and 3 x 10 18 cm' 2 to produce the absolute 0 2 line intensities measured. 0 2 column densities 
this high are hard to produce in quiescent gas (e.g.,|Hollenbach et al.;j 2009j) . Second, Peak 1 is 
the site of intense shock, rather than quiescent cloud, emission. Third, if a component of quies¬ 
cent gas close to H 2 Peak 1 were responsible for the 0 2 emission, and the 0 2 abundance were 
close to its quiescent cloud peak of ~ 5 x 10“ 8 , as in Rho Ophiuchi, the required H 2 column 
density, iV(H 2 ) = iV(0 2 ) / y (0 2 ), would be on the order of 10 25 cm' 2 , which is much higher than 
observed. Thus, process(es) different from those invoked to explain the 0 2 emission toward Rho 
Ophiuchi are needed to explain the detections toward Orion H 2 Peak 1. 


In addition to the presence of shock activity, Orion is the site of O and B stars that produce 
strong FUV radiation. The FUV field near the Trapezium has been estimated to be G 0 > 10 4 
based upon the total radiation from the Trapezium stars - and the O star 6 l Ori C in particular. 
The intensity of this field is corroborated by the strength of the far-infrared [CII] and [O I] fine- 
structure lines mapped toward the Orion molecular ridge, the strength of several near-infrared 
lines whose intensities have been ascribed to recombinations to highly excited states of Cl, 


and the strength of near-infrared NI lines excited by the fluorescence of UV lines (Herrmann 
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et al.||1997} |Marconi et al.||1998| |Walmsley et al.||20 00). The amount by which the FUV field is 


attenuated between the Trapezium and H 2 Peak 1 due to intervening material is uncertain, as 
is the FUV radiation from other B stars in the Orion cloud. However, as discussed in §3, even 
modest amounts of FUV radiation (i.e., G 0 < 10) can affect both the structure of and chemistry 
behind (7-type shocks and motivate the study here. 


FUV-illuminated C-type shocks provide a natural explanation for the Orion H 2 Peak 1 0 2 emis¬ 
sion for three reasons. First, as discussed in §3.2 and §3.3, the FUV field can increase the 
atomic oxygen abundance in the preshock gas, which will increase the peak H 2 0 abundance in 
the postshock gas. Second, the postshock OH abundance (and, ultimately, the 0 2 abundance) is 
increased via the photodissociation of H 2 0 as well as OH-producing chemical reactions, some 
with high activation barriers, enabled by the elevated postshock gas temperatures. By raising 
the 0 2 abundance above that attainable in cold quiescent gas, the implied H 2 column density can 
be brought closer to observed values. Finally, the similarity between the LSR velocities mea¬ 
sured for the 0 2 lines of ~ 11 km s -1 (see Chen et al.), and those measured for the 22 GHz and 
621 GHz H 2 0 masers of 10- 13 km s -1 , suggests a possible physical connection. As discussed 
in §4.1, the maser emission indicates that some shocks are propagating in the plane of the sky, 
which can potentially provide a higher line-of-sight column density of 0 2 than face-on shocks 
while allowing for the narrow observed 0 2 line widths. 


For these reasons, as well as the conspicuous presence of shock activity associated with H 2 Peak 


1, both [ChenetaL (2014) and we focus on shocks as the most likely source of the detected 0 2 
emission. However, three significant differences distinguish the approach previously taken by 
Chen et al. and that taken here. First, Chen et al. assume that the preshock gas and the shock are 
fully shielded from FUV radiation, with only the postshock gas subject to a FUV field. Here, 
we assume that the preshock, shock, and postshock gas are illuminated by a common FUV field. 
Given that the width of the shock front is approximately 10 16 cm or less (see Figs. [9] [10| and 
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SD- this approach seems well justified. In particular, the analysis presented here avoids the 
inconsistency in which the postshock physical conditions are determined assuming the absence 
of FUV radiation while the postshock chemistry, which is governed by these physical conditions 
(e.g., density and temperature), requires the presence of FUV radiation. 


Second, the above is important since the effects of FUV radiation on the shock structure can 
result in substantial changes to the shock width and postshock temperature, with the latter af¬ 
fecting the postshock chemistry. As shown in Fig. [2j for a given preshock density and shock 
velocity, increasing the FUV field intensity both increases the peak postshock temperature and 


reduces the shock width. As shown in Bergin etal. (1998), the formation rate of H 2 0 via the 
neutral-neutral reactions 0+H 2 —> OH+H; OH+H 2 —» H 2 0+H becomes important when the 
postshock temperature exceeds ~ 300 K and increases rapidly with temperature. In fact, as the 
peak postshock temperature rises with increasing G 0 , the H 2 0 formation rate begins to exceed 
the H 2 0 photodestruction rate resulting in a net increase in the postshock H 2 0 abundance with 
G 0 (e.g., compare Fig. [2] middle and bottom panels). This trend continues up to the full con¬ 
version into H 2 0 of all gas-phase O nuclei not locked in CO. Only when the postshock cooling 
causes the temperature to drop below about 300 K does the H 2 0 formation rate decrease signifi¬ 
cantly. When this occurs, photodestruction dominates formation and the H 2 0 abundance begins 
to decline, as shown in Fig. |2j Ignoring the effects of FUV radiation leads to an underestimate 
of the postshock temperature and the postshock depth to which H 2 0 exists in high abundance 
as well as an overestimate of the shock width, all of which affect the resulting 0 2 emission. 


Third, since 0 2 is formed from the photo-destroyed products of H 2 0, i.e., O and OH, H 2 0 must 
be relatively abundant (i.e., \ (H 2 0) > 10” 5 ) in the postshock gas to produce the observed 0 2 
emission. There are two ways to achieve this: (1) by placing H 2 0 into the gas phase directly 
through the sputtering of H 2 0 from the ice-covered mantles of dust grains behind the shock; 
and, (2) by chemically forming H 2 0 behind the shock via the neutral-neutral reactions above. 
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However, for process (2) to form a high H 2 0 abundance, the preshock gas must have a high 
abundance of atomic oxygen. Within gas in which the H 2 0 depletion onto grains is significant, 
Chen etal. require that the sputtering efficiency be > 10% behind a 12 km s 1 shock. This 


efficiency exceeds all theoretical predictions of which we are aware, including Draine etal. 


( | 1 983[ ), [Jlmenez-Serra etal.| ( f2008[ see their Figure 7), and |Flower & Pineau des Florets | ( ]20 1 0] ) . 


Recently, Neufeld etal. (2014) have modeled the water abundance behind interstellar shocks 
based on Herschel measurements of far-infrared (IR) and submillimeter measurements of CO 
and H 2 0 in combination with Spitzer measurements of mid-IR H 2 rotational emission. Their 
best-fit results are in good agreement with the prediction that only when shocks reach a velocity 
of > 25 km s' 1 will H 2 0 be completely removed from grain mantles. Thus, it is not clear that 
sufficient H 2 0 can be removed from grain mantles to support the model offered in Chen et al. if 
depletion of H 2 0 onto dust grains is significant in the preshock gas and negligible sputtering 
efficiency is assumed behind a 12 km s" 1 shock. 


Alternately, Chen et al. suggest that the preshock gas is sufficiently young that depletion is min¬ 
imal, in which case the atomic oxygen abundance is high. However, since this model assumes 
the preshock gas is collapsing from an initial density of 100 cm~ 3 , the time between when the 
density evolves to its required preshock value (i.e., 4.2 x 10 4 cm~ 3 ) and when depletion be¬ 
comes important is ~ 1.7 x 10 6 years. Nevertheless, invoking a short interval during which the 
preshock is sufficiently dense, but undepleted, could result in the high preshock O abundance 
needed to yield a postshock H 2 0 abundance sufficient to produce detectable 0 2 emission. A 
short-lived period in which both the preshock gas density and O abundance are high could also 
account for the rarity of 0 2 detections, as discussed by Chen et al. However, this model still 
presumes the preshock gas is fully shielded while the postshock gas is not, an assumption not 
justified in the Chen et al. model. 


Introducing a FUV field to the preshock gas solves these problems. FUV photons can photodes- 
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orb H 2 0 from grain mantles when shock velocities are too low to produce significant sputtering 
(i.e., v s < 20 km s _1 ). (When shock velocities are > 25 km s _1 , sputtering will remove any H 2 0 
not photodesorbed in the preshock gas.) Even if the FUV field results in the photodestruction of 
gas-phase H 2 0 in the preshock gas, the O and OH produced will be rapidly converted to H 2 0 in 
the warm postshock gas. Thus, a FUV field can eliminate the need for high assumed sputtering 
efficiencies at low shock velocities while relaxing the need to invoke a particular epoch in the 
time-dependent evolution of the preshock gas. As for time-dependent preshock evolution being 
the cause for the rarity of 0 2 detections, a different explanation is suggested below. 

Table 2 lists shock models considered here that both reproduce the absolute 0 2 line intensities 
and line ratios for incident FUV fields having a G 0 = 0.1, 1, and 10. Several trends are evident. 
First, as discussed in §3, as G a increases, the ion-neutral coupling increases and the width of 
the shock decreases. This is evident in the scale of the x-axcs in Figs. [2} [9j [T0| |TT] Fikewise, 
as G 0 increases, roughly the same postshock temperatures are achieved with decreasing shock 
velocities. Second, as G 0 increases, the increased photodissociation of molecules lowers the 
peak postshock abundance of H 2 0, OH, and 0 2 , thus requiring greater H 2 column densities to 
produce the same optically-thin 0 2 line flux (see Table 2). 


Within the observational uncertainties, all three models in Table 2 can explain both the measured 
0 2 line fluxes and line ratios; is one model preferred? One possible discriminator is the degree 
to which each model agrees with estimates of the total H 2 column density and number density 


toward Peak 1. Rosenthal, Bertoldi & Drapatz (2000) estimate the H 2 column density to be « 
2 x 10 22 cm -2 based on the observation of 56 H 2 ro-vibrational and pure rotational lines toward 


Peak 1 using the Infrared Space Obsen’atory’s Short Wavelength Spectrometer (de Graauw 


etal.||1996 ). However, because the lowest-lying H 2 pure rotational transition lies 510 K above 
the ground state, estimates based on H 2 lines are more sensitive to the warm and hot gas and 
much less so to the cooler (i.e., T < 100 K) gas where the 0 2 emission peaks. 
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The cooler gas is better, if less directly, sampled by 12 CO and its isotopologues. Using multi- 
line transitions of 12 CO, 13 CO, and C 18 0, Peng etal. (2012) derive a 12 CO column density of 


between 5 x 10 18 cm 2 and 10 iJ cm 2 toward Peak 1, in good agreement with earlier esti¬ 


mates of Gonzalez-Alfonso etal. (2002). This corresponds to an H 2 column density of between 


6 x 10 22 cm 2 and 1.2 x 10 23 cm 2 , assuming a 12 CO abundance of 8 x 10 5 ([Wilson & Mat- 


teucci 1992). However, even these studies rely on transitions with upper level temperatures 
above the ground state > 32 K and may exclude contributions from a colder component. 

Estimates of the preshock H 2 number density vary depending on the tracer being modeled. 


Generally, estimates exceed 10 5 cm~ 3 , and can range as high as 10 8 cm~ 3 (cf. Kaufman & 
Neufeld 1996a Rosenthal, Bertoldi & Drapatz |2000[ Gonzalez-Alfonso etal. 2002). None of 
the models we are aware of suggest that preshock densities toward Peak 1 are less than 10 4 cm' 3 . 


While none of the models we considered is in perfect agreement with the Peak 1 H 2 column 
density and number density estimates, the best-fit model in Table 2 that assumes G 0 — 1 comes 
closest. Moreover, the shock velocity of this model (i.e., 23 km s _1 ) is within the range of 


typical shock speeds in OMC-1 of 30± 10 km s 1 (Kristensen etal. 2003 2007). 


The best-fit G 0 = 0.1 model is a good approximation to the fully-shielded (7-type shock models 
that have traditionally been used to model a wide variety of molecular emission toward outflows. 
We show that fully-shielded (i.e., G 0 < 0.1) models fail to reproduce the 0 2 line intensities and 
ratios with shock conditions inferred from other data; however, allowing for the presence of 
even modest amounts of FUV radiation (i.e., G 0 > 1) provides a better fit to the observations. 


Finally, we note that the circumstances that explain the 0 2 emission toward Orion H 2 Peak 1 
may also explain the rarity of similar detections in other sources. First, the shocked gas emission 
around Peak 1 subtends ~ 20", larger than for most outflows. Since the strength of the optically - 
thin line emission within the beam scales with the 0 2 column density x the area of the emitting 
region, the column density requirements are eased if the emitting area is large, as toward Peak 
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1. Second, the presence of strong H 2 0 maser emission toward the position of the 0 2 detections 
is evidence not only for the presence of shocks, but for some shocks propagating in the plane 
of the sky. High postshock temperatures promote the production of 0 2 well above abundance 
levels predicted for cold quiescent gas, while the geometry (in the plane of the sky) provides 
the observer with a potentially large line-of-sight column of gas with this relatively high 0 2 
abundance. Of course, shocks of the sort described here would be expected to produce enhanced 
optically thin 0 2 emission regardless of whether these shocks were propagating perpendicular 
to our line of sight or not. However, if such shocks produce 0 2 integrated intensities comparable 
to what was detected toward Peak 1, but with propagation angles closer to our line of sight, the 
0 2 line profiles would appear broader. Thus, the same total line flux would be spread over a 
greater number of velocity bins with a corresponding reduction in the line center amplitude. For 
weak emission, such as that exhibited by the observed 0 2 lines, narrow lines are generally much 
easier to detect than lower amplitude broad lines, particularly in the presence of noise. The rarity 
of 0 2 detections may therefore result from the need to simultaneously satisfy four conditions: 
(1) an 0 2 emitting area that fills a non-negligible fraction of a beam; (2) the presence of C- 
type shocks having peak postshock temperatures sufficient to drive the efficient production of 
H 2 0; (3) illumination by FUV radiation that can both enhance the preshock gas-phase O-nuclei 
abundance and photodissociate postshock H 2 0; and, (4) a shock propagation vector close to 
perpendicular to the line of sight that naturally provides narrow line profiles. 


We wish to acknowledge helpful discussions with Paul Goldsmith. This work was supported by 
NASA through Astrophysics Data Analysis Program grant No. NNX13AF16G. 
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TABLE 1. Summary of O 2 Observations Toward Orion FL Peak 1 


Rest 

Frequency 

(GHz) 

Transition 

(N, 

Energy Above 
Ground State 

( E u /k ) 

Integrated 0 
Intensity 
(K km s -1 ) 

Line“ 

Width 

(km s" 1 ) 

LSR 
Velocity 
(km s -1 ) 

HIFI Beam 6 
Diameter 

(arcsec) 

487.249 

3, 3 - 1,2 

26 K 

0.081 

3.1 

10.2 

44.7 

773.840 

5,4 - 3,4 

61 K 

0.154 

1.7 

11.0 

28.2 

1120.715 

7, 6 - 5, 6 

115 K 

0.043 

0.8 

11.0 

18.9 


After 


Chen etal. 


(2014); b See 


Roelfsema et al. 


( 2012 ). 
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Log 10 H 2 Density (cm' 3 ) 


Fig. 1. Preshock gas-phase abundances of atomic oxygen (red curves) and H 2 0 (blue curves) as 
a function of preshock H 2 density and FUV field strength, G 0 . The abundances are determined 
from the PDR model used here (see text) assuming a visual extinction A v =1.0 between the FUV 
source and the preshock gas. For a given density, lower values of G 0 result in greater amounts 
of O incorporated into ice mantles and thus not available in the gas phase. In the absence of 
sputtering, the preshock gas phase O abundance sets the maximum postshock H 2 0 abundance. 
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Distance Behind Shock Front (cm)/10 15 Distance Behind Shock Front (cm)/10 15 


Fig. 2. Profiles of postshock gas temperature (left panels) and resulting H 2 0, OH, and 0 2 
abundance (right panels) for G 0 = 0.1 (top row), G 0 = 1 (middle row), and G 0 = 10 (bottom 
row). All models have been run for a preshock H 2 density of 10 5 cm” 3 and a shock velocity 
of 10 km s _1 . Because this shock velocity is low, the preshock gas-phase O abundance is not 
increased by sputtering. The fiducial temperature of 1000 K is denoted by a horizontal dashed 
line in the left panels to highlight the increase in peak postshock temperature with increasing G a . 
Also note that distance behind the shock front at which the gas temperature and H 2 0, OH, and 
0 2 abundances peak, and the length scale over which the elevated temperatures and abundances 
persist, decrease with increasing G a . 




















Page 32 



(ooosr) 'oeo 


Fig. 3. Position of 0 2 487 GHz ( red circle), 11 A GHz ( green circle ), and 1121 GHz ( blue circle ) 
beams centered on H 2 Peak 1 (a 2 ooo = 5 h 35 m 13?7, 5 2 ooo = —5° 22' 09"; after Chen etal. 2014) 
superposed on the H 2 image obtained by Cunningham ( 2006| ). The Herschel/HWl instrument 
was also used to map a ~ 2' x 2 ' region centered on BN/KL; only those positions toward Peak 1 
exhibit evidence for H 2 0 maser emission ( |Neufeld et al.||2013| ). The 621 GHz HIFI beams and 
their sky positions are shown as filled yellow circles. The H 2 0 spectra, with the prominent 
maser feature at « 11 km s -1 , is shown on the right. 
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Fig. 4. Proposed scenario in which the 0 2 emission arises within gas in which the 0 2 abundance 
is enriched by the passage of a non-dissociative (7-type shock. If these shocks are propagating 
perpendicular to our line of sight, as suggested by the presence of H 2 0 masers toward Orion 
H 2 Peak 1, the column density of warm, 0 2 -enriched gas can be significantly greater than for 
face-on shocks. Such a scenario permits a line-of-sight LSR velocity shift of ~ 12 km s' 1 , 
while the perpendicular shock velocity can be a free parameter. Finally, because the line of 
sight samples gas with a low velocity gradient, consistent with maser production, the measured 
0 2 line widths of <3 km s _1 can be understood. The models described here divide the postshock 
region into 1,000 zones, each of which is characterized by a distinct density, temperature, and 
set of chemical abundances. The sum of the computed 0 2 emission from all of these zones is 
then compared with the observed integrated intensity in each of the 0 2 lines detected. 
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Fig. 5. Plot of the agreement between the predicted and observed 487 GHz, 774 GHz, and 
1121 GHz 0 2 line ratios produced by various shock models versus H 2 column density and 
shock velocity for G 0 = 0.1. As discussed in the text, the lower the value of the Goodness of 
Fit, the better the shock conditions reproduce the observed 0 2 line ratios, with a value of zero 
representing a perfect fit. Logi 0 of the preshock H 2 density is shown for each family of models. 
The required H 2 (and 0 2 ) column density scales inversely with the area of the 0 2 emitting 
region. The above plot assumes a 487 GHz 0 2 emitting area of 400 sq. arcseconds (see text). 
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Fig. 6. Plot of the agreement between the predicted and observed 487 GHz, 774 GHz, and 1121 
GHz 0 2 line ratios produced by various shock models versus H 2 column density and shock 
velocity for G 0 = 1. As discussed in the text, the lower the value of the Goodness of Fit, the better 
the shock conditions reproduce the observed 0 2 line ratios, with a value of zero representing a 
perfect fit. Logi 0 of the preshock H 2 density is shown for each family of models. The required 
H 2 (and 0 2 ) column density scales inversely with the area of the 0 2 emitting region. The above 
plot assumes a 487 GHz 0 2 emitting area of 400 sq. arcseconds (see text). 
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Fig. 7. Plot of the agreement between the predicted and observed 487 GHz, 774 GHz, and 
1121 GHz O 2 line ratios produced by various shock models versus H 2 column density and 
shock velocity for G 0 = 10. As discussed in the text, the lower the value of the Goodness of 
Fit, the better the shock conditions reproduce the observed 0 2 line ratios, with a value of zero 
representing a perfect fit. Log 10 of the preshock H 2 density is shown for each family of models. 
The required H 2 (and 0 2 ) column density scales inversely with the area of the 0 2 emitting 
region. The above plot assumes a 487 GHz 0 2 emitting area of 400 sq. arcseconds (see text). 
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Fig. 8. Ratio of the observed O 2 integrated line intensities to that of the 487 GHz line (black 
square) given by Chen etal. ( 2014[ ), and those predicted by the best-fit shock model with 
G 0 = 0.1 (red circle), G 0 = 1 {green triangle), and G 0 = 10 {purple diamond) and detailed 
in Table 2. As discussed in the text, the shock models have been computed to reproduce the 0 2 
487 GHz integrated intensity. The 1 a uncertainties in ratios of the 774 GHz and 1121 GHz line 
intensities to the 487 GHz line intensity are shown. 
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Distance Behind Shock Front /10 16 cm 


Fig. 9. Profiles of integrated intensity {top), abundance {middle), and gas temperature {bottom) 
as a function of distance behind the shock front for the best-fit shock parameters for G 0 = 0.1. 
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Distance Behind Shock Front /10 15 cm 


Fig. 10. Profiles of integrated intensity (top), abundance (middle), and gas temperature (bottom) 
as a function of distance behind the shock front for the best-fit shock parameters for G Q = 1. 
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Distance Behind Shock Front /10 15 cm 


Fig. 11. Profiles of integrated intensity {top), abundance {middle), and gas temperature {bottom) 
as a function of distance behind the shock front for the best-fit shock parameters for G 0 = 10. 


























